Non-equilibrium dynamic critical scaling of the quantum Ising chain 
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We solve for the time-dependent finite-size scaling functions of the ID transverse- field Ising chain 
during a linear- in- time ramp of the field through the quantum critical point. We then simulate 
Mott-insulating bosons in a tilted potential, an experimentally-studied system in the same equilib- 
rium universality class, and demonstrate that universality holds for the dynamics as well. We find 
qualitatively athermal features of the scaling functions, such as negative spin correlations, and show 
that they should be robustly observable within present cold atom experiments. 
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The Kibble-Zurek (KZ) mechanism describes the dy- 
namics of a system as it is ramped across a phase tran- 
sition at a finite rate. Kibble first introduced this idea 
to model symmetry breaking during cooling of the early 
universe [1]. Later, Zurek showed that condensed matter 
systems exhibit similar behavior in the context of slowly 
ramping across the superfluid phase transition in liquid 
^He [2]. Both proposals looked at ramps from deep in 
the disordered phase to deep in the ordered phase. His- 
torically, this has been the primary focus of research on 
the KZ mechanism (c.f.e. [3, 4]), although recent work 
has explored ramps that end at the critical point [5-7]. 

However, the scaling theory of the KZ mechanism 
applies more generally than these two limits. Near 
the quantum critical point (QCP), observables measured 
during the ramp are postulated to exhibit universal scal- 
ing forms [8, 9]. In this Letter we solve for the scaling 
functions of excess heat and equal-time spin correlations 
in the ID transverse-field Ising (TFI) chain for two scal- 
ing directions - time t during the ramp and finite size L 
- and investigate the qualitative features of these observ- 
ables. In many cases, the quantum dynamics of a closed 
non-integrable quantum system lead to a thermal state. 
In contrast, we show here that the spin correlation func- 
tions are qualitatively different from those of any thermal 
state. We also find a ramp protocol where the state at 
long time not only does not thermalize, but also does 
not dephase to the diagonal ("generalized Gibbs") en- 
semble. We provide evidence for the universality of the 
dynamics by using time-dependent matrix product states 
(tMPS) to simulate the experimentally-realizable, non- 
integrable model of Mott insulating bosons in a tilted 
potential [10, 11]. We see that the athermal features of 
the TFI chain are robust for small systems and discuss 
the prospects of seeing scaling collapse in present-day ex- 
periments. 

Transverse- field Ising chain - The Hamiltonian for the 
TFI chain on an L-site ID lattice is 

L 



where A is a tunable transverse magnetic field and s^'^ are 
Pauli matrices with periodic boundary conditions (si = 
We work in units where h = 1 and J = 1/2. 
This model has a quantum phase transition (QPT) at 
Ac = [12]. In equilibrium, the system is in a disordered 
paramagnetic (PM) phase for A < 0, while for < A < 2, 
the system is in a ferromagnetic (FM) phase with two 
degenerate ordered ground states in the limit L ^ oo. 

The Hamiltonian of the TFI chain can be diagonalized 
by applying a Jordan- Wigner transformation [12, 13]: 
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where Cj creates a fermion at site j. As the Hamiltonian 
conserves fermion number modulo 2 and the ground state 
has an even number of fermions [13], we restrict ourselves 
to that subspace. Fourier transforming. 
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This Hamiltonian can only excite fermions in momen- 
tum pairs (/c, — /c), so we introduce a pseudo-spin a cor- 
responding to whether the (/c, —k) pair is filled (a^ = 1) 
or unfilled (<j^ = —1). Then in the sector where all 
fermions are paired, which contains the ground state, 
H = Yll=QHk, where 



Hk = (I - \ - cos k)(j^ + (sin k)(j'^ . 
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Note that is a free-fermion Hamiltonian, so the TFI 
chain is integrable. 

In a KZ ramp, A is varied as a function of time near 
the critical point. For simplicity, we focus primarily on 
the case of a linear ramp, \{t) = vt^ where v is the ramp 
rate and the QCP is at A = t = 0. The ramp begins at 
t = — oo deep in the disordered phase, where the wave- 
function is initialized in the ground state of the instanta- 
neous Hamiltonian H{t = — oo). For an infinite system, 
the Hamiltonian is gapless at the critical point. There- 
fore, it is impossible to ramp slowly enough through the 
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QCP to avoid creating excitations and to produce true 
long range order in the ordered phase. 

Near the QCP, critical slowing down tells us that the 
characteristic time scale ~ A~^^ becomes arbitrarily 
large, where v and z are the (positive) correlation length 
and dynamic critical exponents respectively. Thus for 
a non-zero ramp rate, there exists a Kibble time, —Ik^ 
at which the lowest momentum modes stop following 
the ramp adiabatically and become significantly excited. 
More explicitly. 
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For the TFI chain = z = 1, so ~ v One 
can similarly define a length Ik ~ ^//^ such that Ik and 
Ik set the characteristic time and length scales for the 
Kibble-Zurek scaling forms [8, 9]. 

For a KZ ramp of the TFI chain, the wavefunction 
can be written as a product = (8)/c|^/c(^)), 

where each mode evolves independently as 

= [(1 -vt- COS k)crl + (sin k)cr^] \^k) -(6) 

The KZ scaling limit is defined as taking v ^ with r, 
A, and constant [8, 9], where 

r tv^^'^ - t/tK 
A = Lv^^^ - L/£k 
= kv-^^'^ - kl 
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Here n = 0, 1, (L/2) — 1 indexes the modes. Note 
that L ^ oo in the scaling limit, as ^ and A is 
constant. In this KZ scaling limit, eq. 6 simplifies to 
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Each mode ^ ^ij) in (8) has the form of a Landau-Zener 
(LZ) equation, which we solve in terms of parabolic cylin- 
der functions [14]. Note that ^/^(t) is expressed solely in 
terms of the scaled variables r and n. 

Observables - A crude measure of deviation from adi- 
abaticity is the excess heat, defined as 
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where |^^(t)) is the instantaneous ground state of H^{t). 
Defining excited state occupancy p^^^ = 1 — |(^^|^^)| 
and energy = 2v^I^\/t^ + ^ the scaled excess heat 
density is 
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For the equal-time s^-s^ correlation function, we pos- 
tulate a non-equilibrium scaling form, and check for 
scaling collapse. The correlation function is defined as 
G{x) = {sjSj_^^)^ assuming translation invariance. In 
equilibrium, has scaling dimension 1/8, so the scaled 
correlation function is 
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where x = 'x>v^l^ . We break up each site into a pair of Ma- 
jorana fermions [12], such that the correlation function 
is the Pfaffian of a matrix whose elements are pairwise 
expectation values of the Majoranas [15]. We evaluate 
the Pfaffian numerically [16], and find that good scaling 
collapse occurs for small v (see fig. la). 

During the initial part of the ramp (r <C —1), the 
system is very nearly in equilibrium (figs. 2a, lb). Deep 
on the FM side of the ramp (r 1), LZ physics tells us 
that the excitation probability of mode n is [14] 
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In between, at finite positive r, the excitation probability 
for each mode - given by the LZ equation - oscillates as a 
function of scaled time before converging to (12). These 
oscillations show up in the excess heat (fig. 2a) and the 
correlation function (fig. lb,c). 

Finite-size effects enter by opening up a gap at the 
QCP. In the regime of small scaled size A < 1, the finite- 
size gap dominates, the system remains near the ground 
state, and the solution can be understood perturbatively. 
In the limit A ^ oo, on the other hand, finite-size effects 
vanish; we will refer to this as the thermodynamic limit of 
KZ scaling (KZ-TDL). Looking at the correlations right 
at the QCP (fig. 2b), we see the finite-size crossover 
from near-equilibrium power-law correlations for A = 5 
to a non-equilibrium exponential decay of correlations for 
A = 30. Similarly, in fig. 2a we see that the scaled excess 
heat q crosses over from being small at A = 5, where 
the system is nearly adiabatic, to much larger values at 
larger A. 

The TFI chain is integrable, so although a KZ ramp 
through the critical point creates excitations, the result- 
ing excited states differ markedly from equilibrium ther- 
mal states at the same energies. From (12), we see that at 
long times the populations of modes with < (log2)/7r 
are inverted, i.e. at a negative effective temperature. 
This leads to qualitatively athermal physics in the scaled 
correlation function, which goes substantially negative 
by r = 5 over a range of scaled distances, as can be seen 
in figs. 2c and Ic. This is qualitatively different from 
any thermal state, which would have a finite correlation 
length but never negative correlations. Similar behavior 
has been seen outside the scaling regime for both slow 



FIG. 1: (color online) (a) Spin-spin correlation function for the TFI chain for a ramp to the QCP (t = 0), showing scaling 
collapse for a wide range of slow ramp rates. The inset shows correlation functions prior to scaling, (b) Scaled correlation 
length and (c) correlation function as a function of scaled time (b,c) and distance (c) during the ramp. The KZ correlation 
length deviates from equilibrium near tv^/'^ = -2.5 and remains finite at the QCP and beyond. We define the correlation 
length as ^v'^^'^ = ^J g{x)x^dx/ 9{x)dx^ where x — xv'^^'^. All data in this figure are in the KZ thermodynamic limit, 
Lv^^^ > 1 (see text). 



[17] and instantaneous [18] quenches of the TFI chain 
and the XY model [19]. 

Finally, we examine what happens formally as r ^ oo 
while remaining in the scaling limit; in particular, does 
the system dephase? For comparison, if one were to 
stop the ramp at some value Xf and wait a very long 
time, the phase differences between modes would increase 
to the point where the phases are essentially random; 
this process is known as dephasing. Dephasing in inte- 
grable systems is a well-studied problem, and it has been 
shown that, in the long time limit, the observables of the 
time-evolved pure state approach those of the generalized 
Gibbs ensemble (GGE, see ref. [20]), given by removing 
all phase information from each mode. Here we define 
the GGE at time r as the dephased ensemble that one 
would approach upon freezing the current Hamiltonian 
and waiting a long time. 
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where |^^(r)) is the excited state of mode Hamiltonian 
H^. In the limit r ^ oo, the mode Hamiltonians asymp- 
tote to (X —cr^^ so the GGE approaches a fixed value 
with 1^0) ^ I t), 1^1) ^ I ;), and ^ e-'^'. 

To see if r ^ oo leads to dephasing, we consider the 
phase difference A(p between characteristic modes = 
and hi = 1^ since excitations are exponentially suppressed 
for hi^l. Starting from some time ^ 1, after which the 
dynamics is effectively adiabatic, the energy difference 
AE and phase difference Acp are 

AE{r) ' ^ ^ ^ 
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Since Acp ^ oo as rj ^ oo, the phase information be- 
tween modes is lost in the long time limit, so the observ- 
ables approach those of the GGE. 

We note that, for non-linear ramps (say cubic ramps, 
A ~ t^), a similar scaling theory can be written down [8]. 
Then the above argument again holds, except now 



A(p : 



2^3 



const. 



(15) 



as T/ oo, implying that phase information remains, 
and the cubic ramp does not approach the GGE [21]. To 
summarize, both linear and cubic ramps exhibit ather- 
mal behavior, such as negative correlations, which come 
from the inversion of low- momentum modes, but only the 
linear ramp dephases to the GGE in the long time limit. 

Tilted bosons - While the TFI chain is a beautiful the- 
oretical model, it is difficult to realize in the lab. How- 
ever, there have been a number of recent experimental 
advances with other, non-integrable models in the Ising 
universality class [11, 24, 25]. Here we focus on one, 
the Mott insulator in a tilted potential (MITP), realized 
experimentally in ref. [11]. The MITP consists of a one- 
dimensional lattice containing a Mott insulator with one 
boson per site. By adding a sufficiently large potential 
gradient (tilt) J, the system undergoes a QPT from 1 
boson per site to alternately and 2 per site, creating 
dipoles on every other bond (see fig. 2a inset). This 
QCP is in the Ising universality class [10], and can be 
described by an effective spin Hamiltonian 
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where / labels the bonds, S^'^ are Pauli matrices residing 
on the bonds, and 7^ is a projector implementing the 
constraint that no two neighboring bonds are both spin 
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Scaled time, uF^ x/L x/L 



FIG. 2: (color online) Results of tMPS simulations for ramping the experimentally-realizable MITP model (see text), which 
is illustrated in the inset to (a). The data show scaling collapse of MITP (colored dots) to the scaling limit of the TFI chain 
(black lines) as a function of scaled time r = tv-^^"^ (a), scaled system size A = Lv'^^'^ (a,b), and scaled length x/L (b,c) [22]. 
We have checked that data for A = 30 is near the A = co limit for the TFI model (not shown), although finite size effects 
can still be seen in (c). The arrow in the inset to (c) indicates the maximum ramp duration that produces athermal negative 
correlations in a twelve-site bosonic chain (data from Ref. 11, see [23] for details). 



up. The energy scale is set hy u = V2w, where w is the 
hopping energy of the lattice bosons; we work in units 
with u = 1. In analogy with the TFI chain, we define the 
parameters A = 5c — S and v = dtX [22], where Sc = —1.31 
is the location of the QCP [10]. The correlation function 
is given by G{x) = (-1)- ((^|^|+,) - 

We simulate a linear ramp via tMPS, as described in 
ref. [7] . Fig. 2 shows calculated observables of the MITP 
model compared to the Ising scaling forms derived ear- 
lier. We see clear scaling collapse in all three scaling di- 
rections: ramp time t, system size L, and distance x. As 
expected, simulations do not collapse as well for smaller 
system sizes. While the dynamic range of the simulations 
is limited, and the scaling collapse for finite systems is 
imperfect, we consider our results to be strong numerical 
evidence for the postulated universality of KZ scaling for 
systems in the Ising universality class. 

We note that, while scaling collapse for the system 
sizes considered is not perfect at tv^^'^ = 5, the corre- 
lation function goes negative as predicted from the TFI 
scaling function. Therefore, the athermal physics of the 
KZ ramp is robust against small system sizes and the 
breaking of integrability, and is a qualitative feature of 
our model that should be visible experimentally. A ma- 
jor open question is whether such athermal behavior will 
manifest in KZ ramps near non-integrable QCPs, where 
the relatively short-time dynamics of the KZ ramp will be 
in competition with the long-time expectation of eigen- 
state thermalization [26]. 

Finally, we compare our time scales to the those of the 
real experimental system [11]. Using the experimental 
parameters w ^ 10 Hz and U ~ 400 Hz, where U is the 
on-site repulsion, we estimate the ramp rate necessary 
to see athermal negative spin correlation; our predicted 
ramp rates are shown in the inset to Fig. 2c compared 
to experimental data from Ref. 11. Clearly the necessary 



ramp rates are well within the experimental range [23]. 

In conclusion, we have solved for the full Kibble-Zurek 
dynamic scaling forms of the one-dimensional transverse 
field Ising chain at zero temperature. We provided nu- 
merical evidence for the universality of these scaling rela- 
tions via tMPS simulations of the MITP model. Finally, 
we determined that the relevant time scales for seeing 
these effects experimentally are within the reach of cur- 
rent technology. To see full scaling collapse the experi- 
mental system sizes need to be larger, but they should 
already be sufficient to see qualitatively similar athermal 
physics. 
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Supplementary information for "Non-equilibrium 
dynamic critical scaling of the quantum Ising chain" 

For theoretical simplicity, in the main text we have 
addressed the behavior of both the TFI chain and the 
MITP model with periodic boundary conditions. In the 
thermodynamic limit, boundary conditions have no effect 
on the bulk behavior. However, in experimental practice 
[11], the MITP model is realized in a finite (indeed, small) 
system with open boundary conditions. To better under- 
stand the experimental consequences of our work, we now 
address the KZ dynamics of the MITP model with open 
boundary conditions. We solve for the spin-spin correla- 
tion function G{x) for the same criteria as in fig. 2 from 
the main text, with the results show in fig. 3. 

Open boundary conditions (o.b.c.) in the MITP model 
differ from periodic boundary conditions (p.b.c.) only by 
the absence of the projector V on the ends of the chain, 
which in p.b.c. projects out configurations with up spins 
on both the first and last sites. As a result, on the ordered 
side of the phase transition, both the first and last sites 
will favor spin up over spin down due to the Zeeman field 
(see fig. 3d). This results in a large even/odd effect with 
regards to the system size: for odd size systems, there 
is no problem creating the staggered AFM state with up 
spins on each end, but for even size systems this ordering 
becomes frustrated. Note that a similar effect can be 
seen with p.b.c: for L odd the AFM is frustrated. 

An additional complication with o.b.c. is that, com- 
pared to the translationally invariant system with p.b.c, 
there is no unique way to define G{x); due to the loss of 
translation invar iance, 

- {S^){S^+,) ^ {S^S]+,) -{S^f, (17) 

where (S^) = ^ Choosing either side of (17) in 

defining G{x) gives quantitatively different results. Sim- 
ilarly, including sites near the end of the chain in the 
definition of G{x) introduces end effects whereas consid- 
ering only spin correlations between sites in the middle 
third of the chain emphasizes the bulk behavior at the 
expense of throwing away data. The results we show in 
this supplement use the definition 

G{x) ^ [i^^^U^) - i^'f] ' (18) 

Empirically (data not shown), the various definitions of 
G{x) that we described above gave quantitatively differ- 
ent, but qualitatively similar, results, that all approached 
the p.b.c solution - as discussed below - in a similar 
manner. 

The spin-spin correlations resulting from a KZ ramp 
of the MITP with o.b.c are shown in fig. 3, panels a- 
c The even/odd effects discussed above are apparent 
in a slow sweep (fig. 3a), with neither even nor odd 
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FIG. 3: (color online) (a-c) Spin-spin correlation function of the MITP model with open boundary conditions, compared to 
the TFI scaling function with p.b.c. (black lines), (d) z components of the individual spins, 5*1, for the same parameters as (c) 
with L = 97, 49, and 13 (inset) sites, where j is the site index. Arrows indicate the 180° phase shift (AFM domain wall) that 
yields the negative correlations in (c). 



size systems matching the p.b.c. solution. In general, 
matching between o.b.c. and p.b.c. solutions improves in 
the limit of fast sweeps (fig. 3b, c). This is not surprising, 
as finite size effects vanish in the KZ-TDL, so boundary 
conditions should have no effect in that limit. Note that 
fig. 3 shows scaled correlations as a function of scaled 
length x/L, and we see that the solutions match better 
at small distances x <C i^, where boundary effects are not 
important, as opposed to xjL ~ 0.5, where boundary 
effects are non-trivial. 

Most importantly, the athermal negative correlations 
that we saw with p.b.c. exist with o.b.c. (fig. 3c), and 
therefore are robust with respect to both boundary con- 
ditions and small system size. In fact, o.b.c. enhances 
the ability to see negative correlations experimentally. As 
seen in fig. 3d, due to the absence of the end projector, 
the system breaks the Z2 Ising symmetry and selects the 
AFM with higher spin-up probability at the ends. By 
breaking the degeneracy of the ground states, the do- 
main walls associated with negative correlations can be 
seen simply by measuring the expectation value of S^^ for 
each site j, rather than measuring correlation functions. 
As expected from fig. 3c, for constant A these domain 
walls occur at a roughly constant value of 2 jL. This pre- 
diction, confirmed in fig. 3d, implies that the position 
of the domain wall relative to the edge of the sample is 
tunable by the ramp rate v. 



Based on these observations, we estimate the ramp 
rates necessary to see athermal spin correlations at dis- 
tance X equal to half the system size, using the definition 
of G{x) from above. The results, shown in the inset to 
Fig. 2c in the main text for a twelve-site bosonic chain, 
indicate the the time scales necessary to see this phe- 
nomenon are already within experimental reach. How- 
ever, we should note that for the smallest system shown 
in Ref. [11], a six-site bosonic chain, the negative corre- 
lations do not appear due to finite size effects. For any 
even length chain larger than six sites, these athermal 
correlations do robustly appear. 

In conclusion, while open boundary conditions have 
a non-trivial effect on the scaling dynamics of the MITP 
model, these effects are diminished in the limit Lv^l'^ ^ 1 
or X <C Crucially, the qualitatively athermal property 
of negative correlations survives for o.b.c. and small sys- 
tem size L. In fact, the breaking of the ground state de- 
generacy that occurs as a result of boundary conditions 
makes seeing these negative correlations easier experi- 
mentally, as domain walls in measurements of . There- 
fore, we believe that many of the qualitative features of 
the TFI scaling functions should be experimentally acces- 
sible, while the quantitative scaling collapse will require 
larger system size or some mechanism to compensate for 
the open boundary conditions (c.f.e. [28]). 



